Glass transitions in two-dimensional suspensions of colloidal ellipsoids 
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We observed a two-step glass transition in monolayers of colloidal ellipsoids by video microscopy. 
The glass transition in the rotational degree of freedom was at a lower density than that in the 
translational degree of freedom. Between the two transitions, ellipsoids formed an orientational 
glass. Approaching the respective glass transitions, the rotational and translational fastest-moving 
particles in the supercooled liquid moved cooperatively and formed clusters with power-law size 
distributions. The mean cluster sizes diverge in power law as approaching the glass transitions. 
The clusters of translational and rotational fastest-moving ellipsoids formed mainly within pseudo- 
nematic domains, and around the domain boundaries, respectively. 



Colloids are outstanding model systems for glass tran- 
sition studies because the trajectories of individual parti- 
cles are measurable by video microscopy pj. In the past 
two decades, significant experimental effort has been ap- 
plied to studying colloidal glasses consisting of isotropic 
particles [ll-Q , but little to anisotropic particles |6i] . The 
glass transition of anisotropic particles has been stud- 
ied in three dimensions (3D) mainly through simulation 
0, Molecular mode-coupling theory (MMCT) pre- 
dicts that particle anisotropy should lead to new phe- 
nomena in glass transitions [^, and some of these 
have been observed in recent 3D simulations of hard el- 
lipsoids [nl [12. MMCT [13, El suggests that hard el- 
lipsoids with an aspect ratio p > 2.5 in 3D can form 
an orientational glass in which rotational degrees of free- 
dom become glass while the center-of-mass motion re- 
mains ergodic [9]. Such a "liquid glass" in analogy 
to a liquid crystal, has not yet been explored in 3D or 
even 2D experiments. Anisotropic particles should also 
enable exploration of the dynamic heterogeneity in the 
rotational degrees of freedom. Moreover the glass transi- 
tions of monodispersed particles have not yet been stud- 
ied in 2D. It is well known that monodispersed spheres 
can be quenched to a glass in 3D, but hardly in 2D even 
at the fastest accessible quenching rate. Hence bidis- 
persed or highly polydispersed spheres have been used in 
experiments |5l llSl. Il6|. simulations 17] and theory |18i] 



for 2D glasses. In contrast, we found that monodispersed 
ellipsoids of intermediate aspect ratio are excellent glass 
formers in 2D because their shape can effectively frus- 
trate crystallization and nematic order. 

Here we investigate the glass transition in monolay- 
ers of colloidal ellipsoids using video microscopy. We 
measured the translational and rotational relaxation 
times, the non-Gaussian parameter of the distribution 
of displacements, and the clusters of cooperative fastest- 
moving particles. These results consistently showed that 
the glass transitions of rotational and translational mo- 
tions occur in two different area fractions, defining an 
intermediate orientational glass phase. 

The ellipsoids were synthesized by stretching poly- 
methyl methacrylate (PMMA) spheres [13, They 



had a small polydispersity of 5.6% with the semi-long axis 
a = 3.33 fiia and the semi-short axes h — c — 0.56 /im. 
3niM sodium dodecyl sulfate (SDS) was added to sta- 
bilize ellipsoids and the > 3 mM ionic strength in the 
aqueous suspension made ellipsoids moderately hard par- 
ticles. A monolayer of ellipsoids was strongly confined 
between two glass walls [l^l . Light interference measure- 
ments showed that the wall separation varied by only 
^30 nm per 1 mm so the walls could be considered 
as parallel within the field of view. The area fraction 
(j) = nabp, where p is the number density averaged over 
all video frames. Twelve densities were measured in the 
range 0.20 < (j) < 0.81. During the three to six hours 
measurments at each (f), no drift flow or density change 
was observed. The thermal motion of the ellipsoids was 
recorded using a charge- coupled device camera resolv- 
ing 1392x1040 pixels at 1 frame per second (fps) for the 
highest five concentrations and at 3 fps for lower concen- 
trations. The center-of-mass positions and orientations 
of individual ellipsoids were tracked using our image pro- 
cessing algorithm 2l|. The angular resolution was 1° and 
the spatial resolutions were 0.12 pm and 0.04 pm along 
the long and the short axes respectively. More experi- 
mental details are in the Supplemental Material (SM). 

At high densities the ellipsoids spontaneously formed 
small pseudo-nematic domains with branch-like struc- 
tures each involving about 10^ particles, see Fig. SI of 
the SM. The translational relaxation was characterized 
by the self-intermediate scattering function Fs{q,t) = 



N „»q.(x,(t)-x,(0)) 



(EL, 



)/N where Xj(t) is the position of 



ellipsoid j at time t, N is the total number of particles, 
q is the scattering vector and ( ) denotes a time aver- 
age. In Fig. [Ija), we chose qm = 2.3 pm~^ measured 
from the first peak position in the structure factor at 
high density. The rotational relaxation can be charac- 
terized by the n"^ order of the orientational correlation 
function L„(t) = (E^i cosn{6j{t) — 9j{0)))/N where n 
is a positive integer and 9j is the orientation of ellip- 
soid j. Ln{t) decays faster for larger n, and different 
choices of n can yield the same glass transition point, 
n = 4 in Fig. [ijb) was chosen so that L„ (i) can be better 
displayed within our measured time scales. At high </>, 
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FIG. 2: (color online) (a) The non-Gaussian parameters 
of translational displacements along the long axis (q:2(*), 
solid symbols) and the short axis (a^(f), open symbols). 
<?!> = 0.70,0.74,0.77,0.81 as labeled in the figures, (b) The 
non-Gaussian parameters of rotational displacements. 



FIG. 1: (color online) (a) The self-intermediate scattering 
function Fs{q,t) at qm ~ 2.3 /im~^ and (b) the orientational 
correlation L4{t) for different area fractions, (c) The exponent 
P of the fitting function e"'*''^' for the long-time Fs{qm,t) 
and L4{t). (d) The fitted relaxation time r((^) ~ {(j>c — 4')~'^ ■ 
Solid symbols: different choices of q in F{q, t) for the transla- 
tional motion. Open symbols: different choices of n in Lnit) 
for the orientational motion. 



both Fs{qm,,t) and Li{t) develop two-step relaxations, 
the characteristics upon approaching the glass transi- 
tion. The short-time /3-relaxation corresponds to motion 
within cages of neighboring particles, and the long-time 
a-relaxation reflects structural rearrangement involving 
a series of cage breakings. According to mode-coupling 
theory (MCT), the a-relaxation follows e"'*/'^''^. Fig- 
ure [Ijc) shows that f3 decreases with density indicating 
dynamic slowing down upon supercooling [ill, ll^ • 

MCT predicts that the relaxation time T{(j)) diverges 
algebraically approaching the critical point (j)c: T{(f)) ~ 
((t)c - (f)-^ where 7 = l/(2a) -I- 1/(26) Here a and 
h are the exponents in the critical-decay law Fs{q,t) — 
fl^ + hqt^°- and the von Schweidler law Fs{q, t) — f^ — hqt'^ 
at the initial stage of the /3-relaxation and the crossover 
time to the a-relaxation respectively. The fitted a or 
h is almost a constant at different 0, indicating that 
Fs{q,t) can collapse onto a master curve in the appro- 
priate time regime. This demonstrates that Fs{q,t) can 
be separated into a q-dependent and a ^-dependent part 
(23} . Interestingly, Ln(t) can similarly collapse. The fit- 
ted ar = 0.3 ± 0.02 and br = 0.63 ± 0.02 for F,{q„^,t) 
and ae = 0.32 ± 0.02 and be = 0.65 ± 0.02 for Uit) yield 
7T = 2.45±0.05 and 79 = 2.33±0.05 for the translational 
and orientational correlations respectively. These values 
are close to the jt = 2.3 measured for 3D ellipsoids jl2|. 
In Fig. [TJd), T"'^/^ is linear in for different choices of 
q and n. Interestingly, all the scalings show that the 
glass transitions are at (j)^ — 0.72 ± 0.01 for rotational 
motion and (j)^ = 0.79 ± 0.01 for translational motion. 
This indicates three distinct phases: liquid (0 < 0.72), 



an intermediate orientational glass which is liquid-like in 
its translational degrees of freedom but glassy in its ro- 
tational degrees of freedom (0.72 < (f) < 0.79), and the 
glass state for both degrees of freedom {(f) > 0.79). 

Besides the extrapolations in Fig.[ljd), the existence of 
the orientational glass phase was verified from the non- 
Gaussian parameters a2[t) — {Ax'^) /{3{Ax'^)^) — 1 of 
particle displacements Ax during time t [1]. In super- 
cooled liquids, the distribution of Ax is Gaussian at short 
and long times because the motions are diffusive, but it 
becomes non-Gaussian with long tails at the intermedi- 
ate times due to cooperative out-of-cage displacements 
[EflHSl- This behavior is reflected in the peak of a2{t), 
see Fig. [2l As increases, the peak rises and shifts to- 
wards a longer time, indicating growing dynamic hetero- 
geneity on approaching the glass transitions. In contrast, 
the glass phase lacks cooperative out-of-cage motions, so 
a2 (t) exhibits no distinct peak and declines with time [ij . 
Such a sharp change has been regarded as a characteristic 
of a glass transition [ij . Figure [2] clearly shows the glass 
transitions at (p^ — 0.72 ± 0.02 for rotational motion and 
at (f)'^ = 0.79±0.02 for translational motion. In Fig.jU^a), 
aj (t) is always greater than the corresponding a^(i), in- 
dicating that the translational relaxations and coopera- 
tive out-of-cage motions are mainly along the long axes 
of the ellipsoids. 

The two glass transitions can be further conflrmed 
from the spatial distribution of the fastest-moving par- 
ticles which characterizes the structural relaxation and 
dynamic heterogeneity [3. In Fig.d the fastest -moving 
8% of the particles is labeled in colors because the non- 
Gaussian long tail of the distribution of Ax{t*) covers 
about 8% of the population. Here t* corresponds to the 
maximum of a2 [ij. Different choices of t and the per- 
centage yield the similar results. Neighboring fastest- 
moving ellipsoids form clusters and are labeled using the 
same color. Here two ellipsoids are defined as neighbors 
if they overlap after being expanded by 1.5 times and 
their closest distance does not intersect a third parti- 




FIG. 3: (color online) The spatial distributions of the fastest- 
moving 8% of the particles (labeled in colors) in translational 
(a, c, e) and rotational (b, d, f) motions. Ellipsoids in the 
same cluster have the same color, (a, b) The same frame at 
(l> ~ 0.70 (supercooled liquid); (c, d) The same frame a,t (j) — 
0.77 (orientational glass); (e, f) The same frame at </!> = 0.81 
(glass) with ~5500 particles. 



cle. In the supersaturated liquid, most fast particles were 
strongly spatially correlated and formed large extended 
clusters, see Figs. [3] and Fig. S3 in the SM. This demon- 
strates the a-relaxation occurs by cooperative particle 
motion in both the translational and rotational degrees 
of freedom: when one particle moves, another particle 
moves closely following the first. The colloidal glasses, 
in contrast, show no discernible a-relaxation, and the 
fastest particles in /3-relaxation are randomly dispersed 
without forming large clusters as observed in the 3D 
glass transition of colloidal spheres [l| . Figure [3] clearly 
depicts three regimes: both the translational and rota- 
tional fast particles are distributed hetcrogeneously with 
large clusters at ^ < 0.72; the rotational fast particles 
are dispersed homogenously while the translational fast 
particles form large clusters at 0.72 < (j> < 0.79; and 
both types of fast particles are dispersed homogenously 
at > 0.79. 

The spatial distributions of translational and rota- 
tional fast-particle clusters were anticorrelated. Fig- 
ures [2Ia,c) show that most translational fast particles 
belonged to a few large ribbon-like clusters aligned with 
their long axes within the pseudo-nematic domains. In 



FIG. 4: (color online) The probability distribution func- 
tions for the cluster size of (a) translational and (b) rota- 
tional fastest-moving particles. The lines are the best fits 
of P{Nc) ~ N''^. (c) The fitted exponents z^" for rota- 
tional motions and /i^ for translational motions. The ver- 
tical dotted and dashed lines represent the glass transitions 
for rotational and translational motions respectively, (d) The 
weighted mean cluster size (Nc) ~ {(j>c — 4>)~'' where 4>c = 0.71 
and S'l' = 0.79. 



contrast, the clusters of rotational fast particles formed 
branch-like structures extending over several small do- 
mains around the domain boundaries, see Fig.[3l[b). Fast 
rotational ellipsoids moved between domains by coop- 
erative rotational motion. This demonstrates that the 
nematic order within a domain facilitates translational 
relaxation while the orientational disorder near domain 
boundaries promotes rotational relaxation. Fast trans- 
lational particles are responsible for the out-of-cage dif- 
fusion, while fast rotational particles are responsible for 
domain transformations such as splitting, merging and 
rotating. All the phases in Figs. [3ja-f ) contain some iso- 
lated fast translational and rotational particles; they are 
mainly distributed at the domain boundaries with ran- 
dom orientations. 

The cluster sizes of the fast particles, A^c, exhibit 
a power-law distribution P{Nc) ^ as shown in 

Figs. HJa, b). The fitted exponents /i for transla- 
tional and rotational motions change dramatically near 
their respective glass transitions see Fig. IDJc). The 
fi^''^ — 2.0 ± 0.2 for supersaturated liquids is close to 
the /i^ = 2.2 ±0.2 estimated for hard spheres [ij and the 
/i^ = 1.9 ± 0.1 for Lennard- Jones particles in 3D 2Jj, 
while the /i^^^ = 3.2 ±0.1 for glasses is close to the 
/i^ =3.1 estimated for hard spheres in 3D fl*!. Hence 
/i 2± 2.5 might characterize such glass transitions in gen- 
eral. Figure IH^d) shows the weighted mean cluster size 
{Nc) = J2NcPiNc)/T,NcP{Nc) U^Mi at different den- 
sities. Both (N^) and (iVj) diverge on approaching 
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the corresponding <f>c: (Nc) ~ (<^c — with fitted 

r]^ = 0.81 and 77^ = 0.75, indicating growing cooperative 
regions of mobile particles. Similar scaling and 77^ s have 
been observed in a Lennard- Jones system 2J|, but the 
mechanism is not clear. 

We did not observed nematic phase or semetic domains 
found in 3D spherocylinders [25| because 1) the ellipti- 
cal shape facilitates particles changing orientation and 
forming branch- like structures at high densities [2^; 2) 
The 5.6% polydispersity promotes glass formation. 3) 
Long-wavelength fluctuations are stronger in 2D than in 
3D, which can more easily break the long-range order as 
described by Mermin- Wagner theorem. Ellipsoids with 
p ^ 6 appeared to be good glass formers, which can e asily 
preempt any isotropic- nematic (IN) phase transition [27| . 
In contrast, the glass transition can be preempted by 
crystallization for p ~ 1 in 2D, or by an IN transition for 
rods with p > 25 in 3D 

All of the measurements consistently showed that the 
glass transitions for ellipsoids with p — 6 confined be- 
tween two walls are at (f>^ = 0.72 for rotational motion 
and at (f)'^ = 0.79 for translational motion. For longer 
ellipsoids with p = 9 {a = 5.9 /zm, b = c = 0.65 /im), 
0^ = 0.60 ± 0.02 and (/j^ = 0.72 ± 0.02 were observed 
in the two- wall confinement. This suggests that the in- 
termediate regime between (j>^ and i^J increases with the 
aspect ratio, which could be the reason why such an in- 
termediate regime has not been observed in previous 3D 
simulations of ellipsoids with small aspect ratios 12 , 2^ . 
We also observed the two-step glass transitions in mono- 
layers of heavy ellipsoid sediment near one wall, but the 
transitions increased by 3% area fraction because of the 
stronger out-of-plane fluctuations. 

We conclude that colloidal ellipsoids in a 2D system 
exhibit two glass transitions with an intermediate orien- 
tational glass. This behavior has been predicted in 3D 
by MMCT but not studied in 2D before. The two glass 
transitions in the rotational and translational degrees of 
freedom correspond to inter-domain freezing and inner- 
domain freezing respectively. The orientational glass 
regime appears to increase with the aspect ratio. Ap- 
proaching the glass transitions, the structural relaxation 
time and the mean cluster size for cooperative motion 
diverge — typical features of a glass transition 0, [2^ 24 1 . 
Interestingly, the translational and orientational cooper- 
ative motions are anticorrelated in space, which has not 
been predicted in theory or simulation. A similar two- 
step glass transition has been observed in a 3D liquid- 
crystal system and explained as the freezing of the orien- 
tations of the pseudo-nematic domains and the freezing 
of the translational motion within domains [2^. Here 
we directly observed the conjectured pseudo-nematic do- 
mains in ref. [2^. These results at single-particle resolu- 
tion shed new light on the formation of molecular glasses, 
especially at low dimensionality. 
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